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We demonstrate a non-equilibrium phase transition in a dilute thermal atomic gas. The phase 
transition is induced by dipole-dipole interactions between Rydberg atoms which are separated by 
less than the transition wavelength to neighbouring states. In the frequency domain we observe a 
cooperative shift of the Rydberg state which results in intrinsic optical bistability above a critical 
number of excited atoms. In the time domain we observe critical slowing down where the correlation 
length diverges with critical exponent a — —0.53 ±0.10. The atomic emission spectrum in the many- 
body phase provides strong evidence for a superradiant cascade. 



Non-equilibrium systems are found throughout nature, 
for example in ecosystems, financial markets and climate 
[l]. The steady-state of a non-equilibrium system is a 
dynamical equilibrium between driving and dissipative 
processes. As a result, the complex dynamics cannot 
be described within the well-established framework of 
equilibrium thermodynamics. Phase transitions in non- 
equilibrium systems are particularly interesting as they 
are a consequence of macroscopic long-range correlation 
in a system with microscopic interactions [2 . Whilst the 
theoretical analysis of dynamical phase transitions has 
provided much insight [3]-[6] , experimental work is still in 
its infancy. Recent work has detailed the emergence of 
spatial structures in a cold atomic gas [7j and dissipation- 
induced correlation in a cold molecular gas [8]. Here we 
demonstrate a non-equilibrium phase transition [9] in a 
dilute optically-driven thermal gas. 

The long-range correlation in our non-equilibrium sys- 
tem arises due to cooperative dipole-dipole interactions 
between atoms. The level of cooper at ivity can be quanti- 
fied through the cooperativity parameter C = A/'A^/47r^, 
defined by the atomic number density M and the transi- 
tion wavelength A. When the number of atoms per cubic 
transition wavelength is large (C ^ 1), phenomena such 
as superradiance p^OHTS] and level shifts [I3l [14] become 
important. For ground state atoms, the study of cooper- 
ative effects is usually limited to dense media p3), solid- 
state systems [13] and Bose gases confined to cavities 
|15j . Here we are able to access the cooperative regime 
in a dilute gas by performing narrowband laser excita- 
tion to highly-excited Rydberg states (161 HZ] where the 
dipole-dipole interaction [18r,21 and optical non-linearity 
[22-28 are many orders of magnitude larger. The fact 
that the transition wavelength to neighbouring Rydberg 
states is of the order of 1 mm enables us to observe a non- 
equilibrium phase transition on macroscopic, optically- 
resolvable length scales. 

The direct observation of a cooperative non- 
equilibrium phase transition is illustrated schematically 
in Fig. 1. The experiment, shown in Fig. 1(a), consists of 
a thermal vapour of caesium atoms confined in a quartz 



cell with an optical path length of 2 mm [29^ . The atoms 
are driven into the 26p3/2 Rydberg state using three ex- 
citation lasers which co-propagate through the cell 
The probe laser, with wavelength Ap = 852.3 nm and 
Rabi frequency l^p, is resonant between the ground state 
6si/2 and the excited state 6p3/2. 
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FIG. 1. (color online) Observation of a cooperative non- 
equilibrium phase transition, (a) The excitation lasers co- 
propagate through a 2 mm cell of caesium vapour. The first- 
order phase transition is exhibited as (b) an abrupt switch 
in the off-axis fluorescence and (c) intrinsic optical bista- 
bility in the probe laser transmission. At a critical control 
parameter the system discontinuously jumps between single- 
body (green / dark grey) and many-body (orange / light 
grey) behaviour. Experimental parameters: ground state 
density pG — 4.3 x 10^^ cm~^, probe Rabi frequency VLy> — 
27r X 110 MHz, coupling Rabi frequency Qc = 27r x 200 MHz 
and Rydberg Rabi frequency ^Ir = 27r x 330 MHz. The con- 
trol parameter, in this case the Rydberg laser detuning, is 
varied around Ar = 27r x —220 MHz. The dashed white lines 
in (b) indicate the 2 mm optical path length. 
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FIG. 2. (color online) Intrinsic optical bistability in the opti- 
cal response of the ensemble, (a) Experimental and (b) theo- 
retical optical response AT as a function of Rydberg laser de- 
tuning Ar for Rydberg Rabi frequency Qr increasing from (i) 
to (iii). Experimental parameters: ground state density pG = 
4.3 X 10^^ cm~^, probe Rabi frequency Qp = 27r x 100 MHz, 
coupling Rabi frequency Qc = 27r x 220 MHz and Rydberg 
Rabi frequency = 27rx (30,80,160) MHz. Theoretical pa- 
rameters: cooperative shift S = 27rx (0,81,153) MHz/AT and 
Gaussian width a = 27rx (60,100,140) MHz. 



quency stabilised to the |65i/2,^ = 4) ^ |6p3/2,^' = 5) 
transition. The coupling laser, with wavelength Ac = 
1469.9 nm and Rabi frequency Qq^ is resonant between 
the excited states 6p3/2 and 7si/2- Using excited state 
polarisation spectroscopy |31j, the laser is stabilised to 
the |6j93/2,F' = 5) ^ \7si/2,F'' = 4) transition. Finally 
the Rydberg laser with wavelength Ar = 790.3 nm and 
Rabi frequency (^r, is tuned around the resonance be- 
tween the excited-state 7si/2 and the Rydberg state 
26p3/2. 

When the number of Rydberg atoms is small {C <C 1), 
the system exists in a few-body disordered state and the 
induced Rydberg dipoles oscillate with a random phase. 
At a critical Rydberg number density, the entire medium 
undergoes a phase transition to a many-body ordered 
state (C ^ 1) where the dipoles oscillate in-phase. In 
our experiment, the control parameter is the number of 
Rydberg atoms in the ensemble and this is determined 
by the ground state number density pc, and the intensity 
/r and detuning Ar of the third-step Rydberg laser. 

The phase transition is observed as a dramatic change 
in the off-axis fluorescence, as shown in Fig. 1(b), clearly 
visible with the naked eye. This longitudinal phase tran- 
sition along the beam axis has been predicted theoreti- 
cally [32] but not previously observed. The phase tran- 
sition can also be observed as a sudden change in the 
optical transmission AT, as shown in Fig. 1(c). Through 



electron shelving [33], the change in probe laser trans- 
mission provides a direct readout of the Rydberg pop- 
ulation. The conditions for switching between low and 
high levels of Rydberg excitation depend on the history 
of the ensemble giving rise to hysteresis and hence optical 
bistability. 

Bistable phenomena occur due to the presence of feed- 
back in a non-linear system [34]. Optical bistability was 
flrst observed by placing a non-linear medium in a Fabry- 
Perot interferometer where the cavity mirrors provide 
feedback [35 . In the case of intrinsic optical bistability, 
which we observe here, the feedback originates from the 
atoms themselves In the cooperative regime, many 
atoms are present within a transition wavelength and are 
therefore coupled by the same dipolar fleld. The resul- 
tant cooperative shift is excitation dependent and leads 
to an effective frequency renormalisation of the resonance 
condition [37] • At a critical level of excitation this coop- 
erative shift results in intrinsic optical bistability [38 . 

The onset of intrinsic optical bistability is shown as 
a function of Rydberg laser detuning Ar for increasing 




12 14 16 18 
Intensity (W/mm^) 



22 



FIG. 3. (color online) Critical slowing down as the tem- 
poral signature of a phase transition, (a) Continuous Ryd- 
berg laser intensity /r scan showing bistability and hystere- 
sis in the optical response AT. (b) Discrete Rydberg laser 
intensity /r scan showing the divergence of the switching 
time to steady-state r around the critical transition inten- 
sity /R,crit ~ 17.5 W/mm^. The switching time diverges as 
(/r — /R,crit)°' with critical exponent a = —0.53 ±0.10 (stan- 
dard deviation error) shown by the dashed line of best fit. 
Ground state density pG = 4.3 x 10^^ cm~^, probe Rabi fre- 
quency Qp = 27r X 150 MHz, coupling Rabi frequency Qc = 
27T X 340 MHz and Rydberg detuning Ar = 27r x -220 MHz. 
The error bars represent the standard deviation error on the 
determination of the laser intensity and switching time. 
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FIG. 4. (color online) Atomic emission spectra in the single-body and many-body phase. The visible fluorescence spectrum 
is shown for (a) pG — 3.1 x 10^^ cm~^ and (c) pG — 4.3 x 10^^ cm~^. In the single-body phase, the spontaneous emission 
originates from high-lying Rydberg states as illustrated in (b). However, in the many-body cooperative phase, the spontaneous 
emission originates from low- lying Rydberg states, as illustrated in (d), due to a superradiant cascade between high-lying 
Rydberg states. The ionisation limits from 6si/2, 6pi/2 and 6p3/2 are shown by thick red vertical lines. The blue shaded 
regions highlight the absence of spontaneous emission between 26p3/2 and ionisation which would occur due to a blackbody 
or coUisional excitation process. The thin cyan vertical lines indicate the dipole-allowed transitions. Probe Rabi frequency 
Qp = 27r X 110 MHz, couphng Rabi frequency 11c = Stt x 200 MHz and Rydberg Rabi frequency ^Ir = 27r x 330 MHz. 



Rydberg Rabi frequency (1r in Fig. 2(a). The Rabi fre- 
quency provides a measure of the coupling strength to 
the Rydberg state and is determined by the transition 
dipole moment and the laser intensity /r. At low Ryd- 
berg Rabi frequency, and therefore low Rydberg popula- 
tion, there is a symmetric increase in probe transmission 
AT around resonance at Ar = 0. At moderate Rydberg 
Rabi frequency (ii), the lineshape becomes asymmetric 
and shifted to lower frequencies, similar to the coopera- 
tive Lamb shift [14 . Finally with high Rydberg Rabi fre- 
quency (iii), the lineshape becomes bistable, with a sharp 
switch between low and high Rydberg population. This 
bistable behaviour is accompanied by hysteresis which is 
dependent upon the direction from which resonance is 
approached (blue (dark grey) curve from higher frequen- 
cies and red (light grey) curve from lower frequencies). 
Importantly, this bistability is measured in steady-state 
and is not a transient phenomena. As a result, within 
the hysteresis window the system can be placed in either 
a single-body or many-body phase for exactly the same 
experimental parameters. 

A simple theoretical model using an excitation- 
dependent shift is shown in Fig. 2(b). The single- 
body symmetrical response in Fig. 2(b)(i) is modelled 
by a Gaussian curve. The asymmetry, bistability and 
hysteresis in the cooperative many-body response in 
Fig. 2(b)(ii-iii) is reproduced by the renormalisation 
Ar Ar — S AT where S is directly proportional to 
the cooperativity parameter. 



The phase transition can be confirmed through the ob- 
servation of critical slowing down. This temporal sig- 
nature of a phase transition occurs as the system ap- 
proaches a critical point and becomes increasingly slow 
at recovering from perturbations [39l[40]. In Fig. 3(a), 
the Rydberg laser intensity /r is varied continuously and 
the phase transition from single-body to many-body dy- 
namics occurs at critical intensity /R,crit ^ 17.5 W/mm^. 
The temporal response of the ensemble is measured by 
discretely varying the Rydberg laser intensity /r and 
measuring the time r to reach steady-state, as illus- 
trated in the inset of Fig. 3(b). At the critical transi- 
tion, the switching time diverges according to the power 
law r oc (/r — /R,crit)" shown by the fitted dashed line. 
The critical exponent a = —0.53 ±0.10 (standard devia- 
tion error) is consistent with previous work on first-order 
phase transitions and optical bistability [4T| |42]. 

The atomic dynamics across the phase transition can 
be further understood by measuring the visible spectrum 
of the off-axis fluorescence shown in Fig. 1(b). The emis- 
sion spectra for the single- and many-body phases are 
shown in Fig. 4(a) and (c), respectively. In the single- 
body phase, the dominant transitions indicated by (i) 
and (ii) involve decay from high-lying Rydberg states to 
the ground states of the s, p and d series. This behaviour, 
highlighted in Fig. 4(b), is consistent with spontaneous 
emission where such transitions dominate due to the uj^ 
dependence in the Einstein A- coefficient. 

In the many-body cooperative phase the emission spec- 
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trum is dramatically modified. The dominant spon- 
taneous emission transitions (i) and (ii) are no longer 
present. Instead, the spontaneous emission now origi- 
nates from a range of low-lying Rydberg states indicated 
by (iii) and (iv) and highlighted in Fig. 4(d). Further- 
more, the absence of emission in the blue shaded regions 
indicates that atoms have not been promoted to higher- 
lying Rydberg states, as would occur in a collisional pro- 
cess. We can also neglect the effects of thermal black- 
body photons because the average number of photons 
per mode at Rydberg-Rydberg transition frequencies is 
much lower than the average number of excited atoms 
[12^. The transitions at 455 nm and 459 nm occur in 
both phases and correspond to decay to 6si/2 from 7p3/2 
and 7pi/2 respectively. 

The emission spectrum in the many-body phase can 
be understood as a superradiant cascade to lower-lying 
Rydberg states [43] . Evidence for a superradiant cascade 
has also been observed in ultracold atoms [44l [45] . When 
the cooperativity on a particular transition is high, the 
atoms emit collectively and in-phase with one another. 
Within the excitation volume, we estimate a Rydberg 
atom number Nr ^ 1 x 10^ [22] and a single-atom spon- 
taneous lifetime r ^ 500 /is on the 26p3/2 to 26si/2 tran- 
sition. We therefore expect that the superradiant decay 
occurs on a timescale Tguper = 't/^r ^ 500 ps. This is 
much faster than the spontaneous transitions indicated 
by (i) and (ii) where r ^ 20 fis. As the cooperativity pa- 
rameter C (X n^, where n is the principal quantum num- 
ber, the superradiant cascade eventually stops and gives 
rise to the observed spontaneous emission from low-lying 
Rydberg states as indicated in Fig. 4(iii) and (iv). 

In summary, we have demonstrated a cooperative non- 
equilibrium phase transition in a dilute thermal atomic 
gas. The observations that have been discussed raise 
interesting possibilities for future non-local propagation 
experiments which utilise the long range cooperative in- 
teraction [46 . Furthermore, this work could be used to 
perform precision sensing [47 around the critical point 
and to study resonant energy transfer [48 on optically- 
resolvable length scales. In addition, studies of the flu- 
orescence in the vicinity of the phase transition could 
provide further insight into the dynamics of strongly- 
interacting dissipative quantum systems [3, 4 . 
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